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Prerequisites for amplicon pyrosequencing of microbial 
methanol utilizers in the environment 



The commercial availability of next generation sequencing (NGS) technologies facilitated 
the assessment of functional groups of microorganisms in the environment with high 
coverage, resolution, and reproducibility. Soil methylotrophs were among the first 
microorganisms in the environment that were assessed with molecular tools, and 
nowadays, as well with NGS technologies. Studies in the past years re-attracted notice 
to the pivotal role of methylotrophs in global conversions of methanol, which mainly 
originates from plants, and is involved in oxidative reactions and ozone formation in the 
atmosphere. Aerobic methanol utilizers belong to Bacteria, yeasts, Ascomycota, and 
molds. Numerous bacterial methylotrophs are facultatively aerobic, and also contribute 
to anaerobic methanol oxidation in the environment, whereas strict anaerobic methanol 
utilizers belong to methanogens and acetogens. The diversity of enzymes catalyzing 
the initial oxidation of methanol is considerable, and comprises at least five different 
enzyme types in aerobes, and one in strict anaerobes. Only the gene of the large 
subunit of pyrroloquinoline quinone (PQQ)-dependent methanol dehydrogenase (MDH; 
mxaF) has been analyzed by environmental pyrosequencing. To enable a comprehensive 
assessment of methanol utilizers in the environment, new primers targeting genes of the 
PQQ MDH in Methylibium (mdh2), of the nicotinamide adenine dinucleotide-dependent 
MDH {mdh), of the methanol oxidoreductase of Actinobacteria (mdo), of the fungal flavin 
adenine nucleotide-dependent alcohol oxidase (modi, mod2, and homologs), and of 
the gene of the large subunit of the methanolxorrinoid methyltransferases (mtaC) in 
methanogens and acetogens need to be developed. Combined stable isotope probing 
of nucleic acids or proteins with amplicon-based NGS are straightforward approaches 
to reveal insights into functions of certain methylotrophic taxa in the global methanol 
cycle. 
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PREFACE 

The commercial launch of pyrosequencing and later on, fur- 
ther next generation sequencing (NGS) technologies for direct 
sequencing of complex PCR amplicons facilitated the assessment 
of microbial communities by geno- and ribotype composition 
with high coverage and resolution, and with a high number of sam- 
ples (e.g., Pilloni et al., 2012). Methylotrophic bacteria were among 
the first group of microorganisms that were analyzed by molecu- 
lar tools in the environment (Holmes etal, 1995). Based on such 
studies substantially more is known on the pivotal role of microor- 
ganisms for global cycles of one-carbon (CI) compounds, such 
as the greenhouse gas methane. Another quantitatively impor- 
tant atmospheric volatile organic compound (VOC) is methanol. 
Nonetheless, in situ activities, environmental drivers, and dis- 
tribution of methanol utilizers in the environment has scarcely 
been addressed (e.g., McDonald and Murrell, 1997), and just few 
years ago, studies re-attracted notice to activities of methanol 
utilizers in the environment (Delmotte etal., 2009; Dixon etal, 
2011). 



METHYLOTROPHS IMPACT ON GLOBAL 0NE-CARB0N 
COMPOUND CYCLING 

Aerobic methylotrophs occur in terrestrial and aquatic environ- 
ments on the whole planet, and have been detected in aerated 
and flooded soils of wetlands, grasslands, tundra, and deserts, and 
occur in the phyllosphere and rhizosphere of plants, in open ocean 
waters and other marine habitats (Giovannoni et al., 2008; Angel 
and Conrad, 2009; Kolb, 2009a; Wieczorek etal, 2011; Bissett 
etal, 2012; Gupta etal, 2012; He etal, 2012; Knief etal., 2012a; 
Vorholt, 2012) suggesting that their unique physiology that allows 
them to utilize reduced CI compounds as carbon and energy 
source is of global relevance in ecosystems. Methylotrophs ubiq- 
uitously occur in terrestrial ecosystems, i.e., likely, since plants 
produce CI compounds. Growing plants emit methanol (up to 
0.1% of the photosynthetic carbon) and traces of chloromethane 
and methane (Keppler etal., 2005; Keppler etal., 2006), dur- 
ing decay of lignocellulosic plant material methanol is released 
(Galbally and Kirstine, 2002), and plant compounds are eventu- 
ally converted to methane under anoxic conditions (Drake et al., 
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2009). CI compounds are highly volatile and thus, are emitted into 
the atmosphere. Consequently, the two most abundant organic 
compounds in the atmosphere are methane and methanol (Forster 
etal., 2007). The steady-state concentration of methanol in the 
atmosphere (1-10 ppb) is about 1,000-fold lower than that of 
atmospheric methane (1,800 ppb; Galbally and Kirstine, 2002; 
Heikes etal., 2002; Jacob etal., 2005). Whereas, the estimated 
global emission rate of methane (~10 Tmol per year) from ter- 
restrial ecosystems is only twice as high as the global terrestrial 
emission rate of methanol (~5 Tmol per year) indicating that 
methanol is substantially more susceptible to atmospheric chemi- 
cal reactions (Galbally and Kirstine, 2002; Jacob etal., 2005; Kolb, 
2009a). Methanol triggers the formation of tropospheric ozone, 
and has indirectly a threefold higher global warming potential 
on a one-hundred-year basis than carbon dioxide (Forster etal., 
2007). 

Release of methane from terrestrial ecosystems into the atmo- 
sphere is reduced by aerobic methylotrophs (Conrad, 1995). 
Many aerated soils in natural ecosystems are even net sinks 
for atmospheric methane, which is often correlated with the 
predominance of certain genotypes, such as USCa (Dunfield, 
2007; Kolb, 2009b). Methanotrophic methylotrophs have been 
addressed in numerous environmental studies by using gene 
markers and other biomarkers, and are one of the most stud- 
ied functional groups of microorganisms in the environment 
(e.g., Dedysh, 2009; Kolb, 2009b; Dorr etal, 2010; Luke etal, 
2010; Luke and Frenzel, 2011). There are more than 400 publi- 
cations on methanotrophs in ecosystems over the past 25 years 
based on keyword searches in literature databases (Web of Knowl- 
edge, 04.07.2013, http://apps.webofknowledge.com) highlighting 
the interest in understanding the role of methanotrophs in the 
global carbon cycle. 

Non-methanotrophic methylotrophs likely have a similar 
importance for the global methanol cycle, a fact that has recently 
been more thoroughly addressed in the phyllosphere, soil, and 
ocean waters (Delmotte etal, 2009; Kolb, 2009a; Knief etal, 
2010a,b; Dixon etal, 2011, 2013; Vorholt, 2012; Stacheter etal, 
2013). The assessment of methanol-utilizing methylotrophs in 
the environment is less straightforward than the detection of 
methanotrophs, since methanol utilizers have a substantially larger 
diversity than methanotrophs, and the enzymes that catalyze the 
diagnostic reaction, i.e., the oxidation of methanol to formalde- 
hyde, are more diverse than methane monooxygenases making the 
detection of non-methanotrophic methanol utilizers more chal- 
lenging (Chistoserdova etal., 2009; Kolb, 2009a; Stacheter etal., 
2013). 

The role of methanol utilizers in global methanol cycling is 
still scarcely investigated and warrants studies that address the 
response, activity, and distribution of methanol utilizers in terres- 
trial and other environments. Hence, suitable gene targets are 
mandatory to analyze methanol-utilizing microorganisms with 
amplicon pyrosequencing or to detect them in metagenomic, 
transcriptomic, and proteomic datasets based on sequence homol- 
ogy. The review will describe the latest knowledge on microbial 
taxa that are capable of methanol oxidation including those 
organisms that putatively utilize methanol under anoxic condi- 
tions, and will identify gene markers that have been and can be 



employed for analysis of PCR amplicons by high-throughput NGS 
techniques. 

FACULTATIVELY AEROBIC METHANOL UTILIZERS 

Microorganisms that have the capability to utilize methanol with 
molecular oxygen as an electron acceptor belong to various phyla 
of Bacteria, and have been found within yeast, mold fungi, 
and Ascomycota (Table 1; Bystrykh etal, 1988; Trotsenko and 
Bystrykh, 1990; Nakagawa etal, 1996; Nozaki etal, 1996; Silva 
etal., 2009; Sipiczki, 2012). Bacterial methanol utilizers belong 
to Alphaproteobacteria, Gammaproteobacteria, Betaproteobac- 
teria, Flavobacteriia, Bacilli, and Actinobacteria. Yet, methanol 
utilization (MUT) among Archaea only occurs in strict anaer- 
obic methanogens. Generally, it is known that several bacterial 
methylotrophs utilize methanol or other CI compounds for dis- 
similation, but cannot assimilate carbon from CI compounds. 
Strain HTCC2181 is a recent example, which demonstrates this 
strategy of CI compound utilization (Giovannoni etal., 2008; 
Halsey etal., 2012). Over 200 aerobic species of methylotrophic 
Bacteria have been described (Tables 1 and 2; Kolb, 2009a). Most 
of the known isolates are Gram-negative. Thus, it is remark- 
able that a second isolate of the genus Bacillus has recently been 
described, which was not enriched on conventional methylotroph 
media suggesting a largely uncovered diversity of Gram-positive 
methanol utilizers in the environment (Table 2; Ling etal., 
2011). 

It is well established that some facultatively aerobic methanol 
utilizers are capable of growth on CI compounds with nitrate 
as an electron acceptor (Kolb, 2009a). In addition, many more 
methylotrophs that have the ability to use nitrate as an alter- 
native electron acceptor have not yet been tested for anaerobic 
methanol oxidation (Bamforth and Quayle, 1978; Kolb, 2009a); 
recent examples, in which the physiology has been thoroughly 
assessed, are Methyloversatilis universalis FAM5, and Methylotenera 
versatilis (Kalyuzhnaya etal., 2012; Lu etal., 2012; Mustakhi- 
mov etal., 2013). In environments with a high nitrogen input 
(for example by fertilization) and turnover, facultative aerobic 
and nitrate-dependent degradation of methanol likely occurs in 
oxygen-limited zones (Lu etal., 2012). Based on the current knowl- 
edge, these organisms are accessible by the same gene markers as 
described in the following section (Figure 1). 

MARKER GENES OF BACTERIAL METHANOL UTILIZERS 

Amplicon-based analysis of the diversity of methanol utiliz- 
ers can be achieved by deep sequencing of genes that are 
diagnostic for methanol oxidation (Figure 1; Stacheter etal., 
2013). The CI metabolism of bacterial methanol utilizers com- 
prises a series of enzymatic reactions, which partially can- 
not be found in other heterotrophs and are thus, diagnostic 
for methylotrophs. The most characteristic enzymatic step is 
the initial oxidation of methanol to formaldehyde (Figure 1). 
The oxidation of methanol can be catalyzed by at least three 
different enzymes in Bacteria. There is a pyrroloquinoline 
quinone (PQQ) -dependent and a nicotinamide adenine din- 
ucleotide (NAD)-dependent methanol dehydrogenase (MDH; 
Devries etal., 1992; Chistoserdova etal., 2009; Krog etal., 
2013). PQQ MDH occurs in Gram-negative Bacteria, whereas 
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Table 1 | Classes and phyla of Bacteria and fungi that contain 
methanol-utilizing methylotrophs based on previous reviews (Kolb, 
2009a; Gvozdev et al., 2012). 



Class/phylum/order 


Representative species 


Actinobacteria 


Brevibacteriaceae 


Brevibacterium casei 


Micrococcaceae 


Arthrobacter methylotrophus 


Mycobacteriaceae 


Mycobacterium gastri 


Nocardiaceae 


Rhodococcus erythropolis 


Pseudonocardiaceae 


Amycolatopsis methanolica 


Bacilli 


Bacillaceae 


Bacillus methanolicus 


Alphaproteobacteria 


Acetobacteraceae 


Acidomonas methanolica 


Beijerinckiaceae 


Methylocapsa aurea 


Bradyrhizobiaceae 


Afipia felis 


Hyphomicrobiaceae 


Angulomicrobium tetraedrale 


Methylobacteriaceae 


Methylobacterium extorquens 


Methylocystaceae 


Methlyopila jiangsuensis 


Phyllobacteriaceae 


Mesorhizobium loti A 


Rhizobiaceae 


Ensifer fredii A 


Rhodobacteraceae 


Paracoccus alkenifer 


Sphingomonadaceae 


Sphingomonas melonis 


Xanthobacteraceae 


Ancylobacter dichloromethanicus 


Betaproteobacteria 


Comomonadaceae 


Variovorax paradoxus 


Methylophilaceae 


Methylophilus glucosoxydans 


Rhodocyclaceae 


Methyloversatilis universalis 


Burkholderiales 


Methylibium aquaticum 


Gammaproteobacteria 


Enterobacteriaceae 


Klebsiella oxytoca 


Methylococcaceae 


Methylococcus capsulatus 


Piscirickettsiaceae 


Methylophaga thiooxydans 


Vibrionaceae 


Photobacterium indicum 


Classification unclear 


Methylohalomonas lacus 


Ascomycota 




Gliocladium deliquescens 




Paecilomyces variotii 




Trichoderma lignorum 


Yeasts 




Candida boidini (and others) 




Hansenula capsulatus (and others) 




Pichia pastoris (and others) 


Mold fungi 




Paecilomyces variotii 




Penicillium chrysogenum 



A Harbors xoxF-like gene mxaP. Growth on methanol has not been tested. 



the NAD MDH is typical of Gram-positive Bacillus strains 
and is encoded by the gene mdh (Chistoserdova etal, 2009). 
Furthermore, in Gram-positive Actinobacteria {Amycolatopsis 
methanolica, Mycobacterium gastri MB19) a methanokNDMA 
(N,N'-dimethyl-4-nitrosoaniline) oxidoreductase (MDO) has 
been reported (Bystrykh etal., 1993; Vanophem etal., 1993; Park 
etal, 2010). 

The gene mxaF encodes the catalytic subunit of PQQ MDH, 
which is composed of two different subunits (i.e., MxaFI). 
However, there is a distantly related homolog in some methy- 
lotrophic Burkholderiales (Methylibium; i.e., the gene was named 
mdh2), which encodes an alternative PQQ-dependent MDH 
(Kalyuzhnaya etal, 2008). Beyond mdh2, a further homolog of 
mxaF is known, i.e., xoxF, for which a functional role in methanol- 
metabolism is under debate. xoxF-like genes (synonymous to 
mxaF') occur in Bradyrhizobiaceae and other rhizobia, and may 
encode functional MDHs (Fitriyanto etal., 2011). Similar genes 
can be frequently detected in soil microbial communities using 
mxaF-specific primers (Table 1; Stacheter etal, 2013). If rhi- 
zobia that do not possess the classical PQQ MDH (i.e., MxaFI; 
Moosvi etal, 2005), also utilize and grow on methanol has 
not systematically been analyzed yet. Nonetheless, Bradyrhi- 
zobium sp. MAFF2 11645 contains a Ce 3+ -inducible XoxF-like 
MDH (Fitriyanto etal, 2011). The first functional proof of XoxF 
as a MDH was demonstrated in the phototroph Rhodobacter 
sphaeroides (Wilson etal., 2008). Studies on Methylobacterium 
extorquens AMI suggest that XoxFl and XoxF2 are involved 
in the regulation of mxaF (Schmidt etal, 2010; Skovran etal., 
2011). Purified XoxFl of Methylobacterium extorquens AMI has 
highest methanol oxidation activities when cells were grown 
with methanol and 30 (jlM La 3+ . These activities were compa- 
rable to the canonical PQQ MDH MxaFI. The purified XoxF 
enzyme contained La 3+ suggesting that XoxF is important as 
a calcium-independent MDH that uses rare earth elements as 
cofactors (Nakagawa etal, 2012). In recently discovered methan- 
otrophs of the phylum Verrucomicrobia, xoxF is the only 
detectable gene that may code for MDH (Op den Camp etal, 
2009). xoxF also occurs in non-methylotrophic bacteria, in 
which its metabolic function is unresolved (Chistoserdova, 2011). 
Thus, the detection of xoxF by NGS in environmental 
gene surveys or their occurrence in metagenomes, tran- 
scriptomes, and proteomes may be a hint to environmental 
methanol oxidation but need to be carefully evaluated based 
on recent and upcoming results from pure cultures of various 
organisms. 

A comprehensive assessment of the genotypic diversity of aer- 
obic methanol utilizers in the environment seems possible when 
mxaF, xoxF-like, mdh2, mdh, and genes of MDO of Actinobac- 
teria are simultaneously analyzed. However, only mxaF has been 
successfully detected to date and PCR primers suitable for envi- 
ronmental surveys of the other genes have not yet been developed 
(McDonald and Murrell, 1997; Neufeld etal, 2007; Stacheter 
etal., 2013). More studies on the function of xoxF in further 
methylotrophs and addressing the physiological role of xoxF in 
organisms that are currently not known as methylotrophs are 
warranted to improve the ability to interprete methylotrophy 
gene-targeting surveys in the environment. The employment of 
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Table 2 | List of methanol-utilizing methylotophs that are not included in a previous survey (Kolb, 2009a). 



Class/phylum 


Isolation source 


fac 


pH 


Reference 


Actinobacteria 


Micrococcus luteus MM7 


Oral Cavity 


+ 


nn 


Anesti etal. (2005) 


Bacilli 


Bacillus vallismortis JY3A 


Soil 


+ 


N 


Ling etal. (2011) 


Alphaproteobacteria 



Ancylobacter dichloromethanicus DM16T 
Ancylobacter oerskovii NS05 
Ancylobacter polymorphus DSM 2457 
Ancylobacter rudongensis JCM 1167 
Ancylobacter vacuolatus DSM 1277 
Methylobacterium bullatum B3.2 
Methylobacterium cerastii C15 
Methylobacterium gnaphalii AB627071 
Methylobacterium goesingense AY364020 
Methylobacterium gossipiicola Gh-105 
Methylobacterium longum DSM 23933 
Methylobacterium marchantiae DSM 21328 
Methylobacterium oxalidis DSM 24028 
Methylobacterium phyilosphaera BMB27 
Methylocapsa aurea DSM 22158 
Methyloferula stellata AR4 
Methlyopila jiangsuensis DSM 22718 
Starkeya koreensis Jip08 
Starkeya novella I AM 12100 
Betaproteobacteria 
Methylobacillus arboreus VKM B-2590 
Methylobacillus gramineusVKM B-2591 
Methylopila musalis MUSAT 
Methylophilus rhizosphaerae BMB147 
Methylophilus glucosoxydans B 
Methylophilus flavus DSM 23073 
Methylophilus luteus DSM 2949 
Methylotenera versatilis JCM 17579 
Methylovorus menthalis DSM 24715 
Variovorax paradoxus 5KTg 
Gammaproteobacteria 
Methylomonas koyamae MG30 
Methylomonas scandinavica SR5 
Methylomonas paludis MG30 
Methylothermus subterraneus DSM 19750 
Methylophaga lonarensis MPL 
Methylophaga sulfidovorans RB-1 
Methylophaga thiooxydans DSMO10 



Soil 
Soil 
Soil 

Rhizosphere 
Soil 

Phyllosphere 

Phyllosphere 

Phyllosphere 

Rhizosphere 

Phyllosphere 

Phyllosphere 

Phyllosphere 

Phyllosphere 

Phyllosphere 

Soil 

Soil 

Activated sludge 

Phyllosphere 

Phyllosphere 

Phyllosphere 

Phyllosphere 

Banana 

Rhizosphere 

Rhizosphere 

Phyllosphere 

Phyllosphere 

Sediment 

Rhizosphere 

Oral Cavity 

Water of rice paddy 

Groundwater 

Torfmoor 

Aquifer 

Sediment 

Sediment 

Marine waters 



+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 

+ 
+ 
+ 

+ 
+ 
+ 
+ 
+ 

+ 
+ 



N 
N 
N 
N 
N 
N 

Act 

nn 

N 

N 

N 

N 

N 

N 

Act 

Ac 

N 

N 

N 

N 
N 
N 
N 
N 
N 
N 
N 
Al 
nn 

Act 
N 

Act 

Act 

Al 

N 

nn 



Firsova etal. (2009) 
Lang etal. (2008) 
Xin etal. (2006) 
Xin etal. (2004) 
Xin etal. (2006) 
Hoppe etal. (2011) 
Wellneretal. (2012) 
Tani etal. (2012a) 
Idris etal. (2006) 
Madhaiyan etal. (2012) 
Knief etal. (2012a,b) 
Schaueretal. (2011) 
Tani etal. (2012b) 
Madhaiyan etal. (2009a) 
Dunfield etal. (2010) 
Vorobevetal. (2011) 
Li etal. (2011) 
Im etal. (2006) 
Im etal. (2006) 

Gogleva etal. (2011) 
Gogleva etal. (2011) 
Doronina etal. (2013) 
Madhaiyan etal. (2009b) 
Doronina etal. (2012) 
Gogleva etal. (2010) 
Gogleva etal. (2010) 
Kalyuzhnaya etal. (2012) 
Doronina etal. (2011) 
Anesti etal. (2005) 

Ogiso etal. (2012) 
Kalyuzhnaya etal. (1999) 
Danilova etal. (2013) 
Hirayama etal. (2011) 
Antony etal. (2012b) 
de Zwart etal. (1996) 
Boden etal. (2010) 



fac, facultative methylotrophic; +, growth in substrates with carbon-carbon bonds; —, no growth in substrates with carbon-carbon bonds; nn, no information in 
species description; enrichment at pH 7; N, neutrophile (growth optimum between pH 6.5 and 8.0); Act, acid tolerant (growth optimum between pH 3.0 und 6.5); Ac, 
acidiphilic (growth optimum below pH 3); Al, alcaliphilic (growth optimum over pH 8). 
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A - Aerobic & Facultative Aerobic Methanol Utilizers 

Proteobacteria 
Verrueomicrobia 



Burkholderiales 




Bacilli 



B- Strict Anaerobic Methanol Utilizers 



Methanogens 



I _ Methanol: com 
noid Methyl 
transferase 



4CH 3 -CoM 



Acetogens 



Methanohcorri 
noid Methyl- 
transferase 




3 CH 3 CO-S-CoA-» 



FIGURE 1 | Model of the CI metabolism of aerobic and facultatively 
aerobic methanol utilizers (A), and of methanol-utilizing methanogens 
and acetogens (B). Various Bacteria that employ PQQ MDH and PQQ 
MDH2 utilize nitrate besides molecular oxygen, and are, thus, facultative 
aerobes. Red, enzymatic reactions indicative for methanol utilization. Gray, 
metabolic crossing points to anabolic pathways. Structural genes are 
targets for the molecular assessment of methanol-utilizing microorganisms 
in the environment. PQQ MDH, PQQ-dependent methanol dehydrogenase; 
PQQ MDH2, alternative PQQ-dependent methanol dehydrogenase in 
Methylibium andBurkholderia strains; FAD AO, FAD-dependent alcohol 
oxidoreductase; MDO, methanol oxidoreductase; NAD MDH, 
NAD-dependent methanol dehydrogenase; FAE1, 
tetrahydromethanopterin-dependent formaldehyde activating enzyme; 
MCH, tetrahydromethanopterin-dependent methenyl-methylene 
cyclohydrolase; FaDH, formaldehyde dehydrogenase; GSH, 
glutathione-dependent formaldehyde oxidation; MySH, 
mycothiol-dependent oxidation of formaldehyde in yeast; CH2 = HF4, 
methylene tetrahydrofolate; CH3-THF methyl-tetrahydrofolate; CH3-C0M, 
methyl-coenzyme M; CH3OH, methanol; CH2O, formaldehyde; COOH 
formate; CO2, carbon dioxide; CO, carbon monoxide; CH3CO-S-C0A, 
acetyl coenzyme A; CH3COOH, acetic acid. Biomass, assimilation of 
carbon occurs at the level of formaldehyde and/or on carbon dioxide in 
aerobic methylotrophs; pathways involved are the Calvin-Benson-Bassham 
cycle, ribulose-monophosphate pathway, and the serine cycle. Assimilation 
of carbon in methanogens is mediated by a unique reductive 
acetyl-CoA-pathway, whereas acetogens form acetyl-CoA as an 
intermediate that can be used for biosnthesis. This figure is based on 
previous articles (Thauer, 1998; Ding etal., 2002; Hagemeier etal., 2006; 
Das etal., 2007; Drake etal., 2008; Chistoserdova etal., 2009; 
Chistoserdova, 2011; Gvozdev etal., 2012). 



genes of MDHs of Gram-positive methylotrophs will enhance the 
environmental detectability of methanol utilizers and will aid to 
understand the role of these largely overlooked methylotrophs for 
methanol conversion in ecosystems. 

In addition to gene markers that are diagnostic for 
methanol oxidation, the genes mch (methenyl:methylene 
tetrahydromethanopterin cyclohydrolase) and/ael (formaldehyde- 



activating enzyme) that are indicative for dissimilatory oxidation 
of formaldehyde by tetrahydromethanopterin-dependent reac- 
tions have been used to detect methylotrophs in the environment 
(Kalyuzhnaya et al., 2004; Kalyuzhnaya and Chistoserdova, 2005; 
Stacheter etal., 2013; Table 3). However, their detection alone 
does not allow for the conclusion that the detected organ- 
isms are methanol utilizers since fael and mch also occur in 
non-methylotrophic heterotrophs (Chistoserdova, 201 1; Stacheter 
etal, 2013). 

ROUTES TOWARD ENVIRONMENTAL DETECTION OF 
METHYL0TR0PHIC YEASTS, MOLDS, AND ASC0MYC0TA 

Fungi employ a unique pathway for methanol oxidation, in 
which methanol is oxidized via formaldehyde and formate to 
carbon dioxide, i.e., the MUT pathway (Hartner and Glieder, 
2006; Figure 1). The initial oxidation of methanol is mediated 
by a flavin adenine nucleotide (FAD) dependent alcohol oxidase 
(FAD AO) that produces formaldehyde and hydrogen peroxide 
(Hartner and Glieder, 2006). FAD AO occurs in various genera 
of yeasts, such as Candida and Pichia, in molds, and Ascomy- 
cota (Table 1; Vandenbosch etal., 1992; Gvozdev etal., 2012). 
Known genes encoding for FAD AOs are modi and mod2, however 
homologs exist of which the function is unresolved (Nakagawa 
etal, 2006; Gvozdev etal., 2012). The use of genes of FAD AO 
for the environmental detection of methylotrophic fungi will 
be still challenging since numerous isoenzymes with likely dif- 
ferent kinetic properties exist (Table 1; Ito etal., 2007). The 
role of the diversity and activity of fungal microorganisms for 
environmental conversion of methanol has scarcely been studied 
(Table 1), and warrants, especially in terrestrial environments, 
future research. 

POTENTIAL GENE MARKERS OF STRICT ANAEROBIC 
METHANOL UTILIZERS 

The quantitative contribution of anaerobic methanol conversion 
in soils has scarcely been analyzed (Conrad and Claus, 2005). 
Beyond facultative aerobic methylotrophs, some strict anaerobes 
utilize methanol, i.e., methylotrophic methanogens and aceto- 
gens. A methanol-utilizing acetogen is Moorella thermoacetica, and 
examples for methanol-utilizing methanogens are Methanosarcina 
acetivorans, Methanolobus sp., and Methanosarcina barkeri (Das 
etal., 2007; Antony etal, 2012a; Penger etal, 2012). Recently, 
the methanol-oxidizing enzyme methanobcorrinoid methyltrans- 
ferase (MtaC), encoded by mtaC, has been structurally charac- 
terized in these organisms (Ding etal, 2002; Hagemeier etal, 
2006). An enzyme with homology to MtaC is upregulated dur- 
ing growth on methanol in methylotrophic acetogens (Zhou 
etal., 2005; Das etal., 2007). Hence, the gene mtaC and its 
homolog in acetogens are promising targets to develop gene- 
based detection of strict anaerobic methanol utilizers in the 
environment. 

ASSESSMENT OF METHANOL UTILIZERS BY AMPUCON 
PYROSEQUENCING 

The advent of NGS technologies allow for a dramatic increase 
of sequence information compared with similar efforts when 
using classic Sanger sequencing (Christen, 2008; Liu etal., 2012). 
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Table 3 | Gene markers of methanol-utilizing microorganisms for amplicon-based pyrosequencing or as targets for homology screens in 
metagenome, -trancriptome, or -proteome datasets. 



Methanol utilizers 


Enzyme 


Function 


Gene marker 


Primers 


Reference 


Pyroseq 0 


Proteobacteria 


POQ MDH 


MeOH ox. 


mxaF 


1003t/1555r 


McDonald and Murrell (1997); 
Neufeld etal. (2007) 


Yes 


Proteobacteria, 


PQQ MDH 


Putatively 


xoxF 


1003f/1555r E 


McDonald and Murrell (1997); 


Yes 


Verrucomicrobia 




MeOH ox. 






Neuteld etal. (2007) 




Burkholderiales 


PQQ MDH2 


MeOH ox. 


mdh2 


Not available A 


Kalyuzhnaya etal. (2008) 


No 


Various Proteobacteria 


FAE 


Other 8 


fael 


taelt/tael r 


Kalyuzhnaya etal. (2004) 


Yes 


non-methylotrophs 
















MCH 


Other 8 


mch 


mch-2a/mch-3 


Vorholt etal. (1999); Kalyuzhnaya 
and Chistoserdova (2005) 


Yes 


Actinobacteria 


MDO 


MeOH ox. 


mdo 


Not available A 


Park etal. (2010) 


No 


Bacilli 


NAD MDH 


MeOH ox. 


mdh 


Not available A 


Devries etal. (1992); Krog etal. 
(2013) 


No 


Methylotrophic 


MethanohCoM 


MeOH ox. 


mtaCf- 


Not available A 


Hagemeieretal. (2006) 


No 


methanogens 


methyl-transferase 
system 












Methylotrophic acetogens MethanohCoM 


MeOH ox. 


mtaCA\ke c 


Not available A 


Das etal. (2007) 


No 




methyl-transferase 














system-like 












Fungi 


FAD AO 


MeOH ox. 


modl, mod2, 
others 


Not available A 


Reid and Fewson (1994); Ozimek 
etal. (2005), Hartnerand Glieder 
(2006); Nakagawa etal. (2006) 


No 



MeOH, methanol oxidation. 

A Primers for this group of genes have not designed and tested in environmental surveys. 

B These enzymes do not oxidize methanol, but are involved in formaldehyde oxidation. These enyzmes also occur in methylotrophs that do not use methanol and in 
non-methylotrophs (Chistoserdova, 2011). 

c Homologs of unknown function are present in methanogens (Ding etal., 2002). 

D Has been used in amplicon pyrosequencing. 

E Detect only mxaF and xoxF-like genes of Proteobacteria. 



Amplicon pyrosequencing (i.e., a synonymous term is pyrotag 
sequencing) is one of the best evaluated and oldest NGS tech- 
nologies (Liu etal, 2012). Long reads of about 700-1000 bp 
are possible, and technically unavoidable sequence errors can 
be removed with established software, such as AmpliconNoise 
(Quince etal., 2011; Rosen etal., 2012). Thus, large datasets 
with over 100,000 reads can be quality filtered, trimmed, 
sorted, and clustered into sequence similarity-defined oper- 
ational taxonomic units (Caporaso etal, 2010; Nebel etal, 

2011) . Amplicon pyrosequencing comes along with higher 
costs per read compared to cheaper technologies, such as 
HiSeq, MiSeq, or Ion Torrent (Liu etal., 2012). However, 
the long-read length of pyrosequencing is especially advan- 
tageous when analyzing amplicons. Beyond that, amplicon- 
based pyrosequencing generates highly reproducible and similar 
community structures when compared to standard community 
fingerprinting techniques, such as terminal restriction frag- 
ment length polymorphism (TRFLP) analysis and thus, can be 
used for reliable genotype composition analyses (Pirloni etal, 

2012) . 



Amplicons can be obtained with primers that contain adapter 
sequences that are needed for emulsion PCR-based amplifica- 
tion to generate nanobead-bound sequence libraries (Ronaghi 
etal., 1998). Usually such primers include a several nucleotide- 
long sequence for the identification of the source of sequence 
(i.e., a barcode), such as an individual sample within the study 
(Pilloni etal., 2012). A consequence is long primers with unspe- 
cific extensions at 5' end. This may lead to reduced sensitivity of 
amplification and to unspecific binding (Berry etal., 2011). One 
strategy to overcome this bias is applying two subsequent PCRs. 
The first PCR is conducted with untagged primers followed by 
a second PCR with tagged primers (Berry etal, 2011). An alter- 
native strategy to minimize these shortcomings is to use primers 
with barcodes, but without adapter sequences. Adapters need to 
be added by ligation after the PCR (Stacheter et al., 2013). A com- 
plication of such an approach is the retrieval of two datasets one 
with sequences starting with forward and one starting with the 
reverse primer, but minimizes bias during amplicon amplifica- 
tion (Stacheter et al., 2013). Environmental detection of mxaF-like 
genotypes of methylotrophs by primers 1003f and 1555r include 
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the detection of xoxF-like genes (Table 3; Stacheter etal, 2013). 
When mxaF '-targeting primers were used, also xoxF-like genes 
were detected in various grassland and forest soils by amplicon 
pyrosequencing (Stacheter etal., 2013) making the interpretation 
of data in regard to the capability of detected microorganisms of 
methanol oxidation more difficult, since the function of xoxF is 
still in part under debate. 

mxaF AND H0M0L0GS FOR ENVIRONMENTAL DETECTION 
OF METHYLOTROPHS 

Analysis of non-methanotrophic methylotrophs by mxaF geno- 
typing has been employed in several studies. Nonetheless, only 
one study exist that employed amplicon pyrosequencing (Table 4) . 
All other amplicon-based NGS studies addressed methanotrophs 
and mostly analyzed pmoA (i.e., encodes a gene of a subunit of 
the particulate methane monooxygenase). The use of amplicon 
pyrosequencing is a great step forward toward complete cover- 
age of the real diversity that exists in a given habitat. In this 
review, the authors argue in favor to target structural genes of 
methanol utilizers. One advantage of the use of structural genes is 
the increased sensitivity since rare groups, such as methylotrophs 
in soil communities, can be more reliably detected than by a 
16S rRNA gene-based survey. Several methylotrophs occur in 
taxa of which only some members are capable of methylotro- 
phy (e.g., Bacillaceae; Tables 1 and 2). The detection of such 
methylotrophs by 16S rRNA genes can be misleading, and thus, 
another advantage of the use of genes encoding a methanol- 
oxidizing enzyme is that the detection of the gene marker is 
linked with the potential phenotype of MUT. Nonetheless, gene 
marker-based phylogenies are not always congruent with organ- 
ismal phylogenies (i.e., due to horizontal gene exchange between 
distantly related bacteria or evolution of functionally slightly dif- 
ferent enzymes in the same organism; Friedrich, 2002). In general, 
mxaF-based phylogenies correlate with organismal phylogenies on 



the level of families of methylotrophs (Kist and Tate, 2013a,b; Lau 
etal., 2013). However, for other genes (mdh2, mdo, mdh, modi, 
mod2, mtaC) of methanol-oxidizing enzymes, congruence with 
organismal phylogenies needs to be evaluated. 

Recent evaluation of phylogenetic resolution of mxaF com- 
pared to organismal phylogenies revealed contradicting results 
(Kist and Tate, 20 1 3a,b; Lau et al, 20 1 3 ) . The congruency with the 
16S rRNA gene phylogeny and the resolution of mxaF is sufficient 
(except for some "anomalies") in the non-methanotrophic genus 
Methylobacterium (affiliates with Alphaproteobacteria), i.e., the 
so-called pink-pigmented, facultatively methylotrophic (PPFM) 
bacteria (Kist and Tate, 2013a,b), and mxaF-based taxonomic res- 
olution might be even higher than that of 1 6S rRNA genes (Kist and 
Tate, 2013b). Nonetheless, strain-level identification is not possi- 
ble and requires the analysis of more variable genomic regions 
(Knief et al., 2008, 2010b). Some alphaproteobacterial genera har- 
bor mxaF-like genes that are similar to those of Methylobacterium 
suggesting the occurrence horizontal gene transfer events in evo- 
lution of methylotrophs; Methylobacterium nodularis ORS 2060A 
carries a plasmid with methylotrophy genes including mxaF, sug- 
gesting that this species has acquired this gene from another PPFM 
bacterium (Kist and Tate, 2013a). Using mxaF as a phylogenetic 
marker of methanotrophic Proteobacteria revealed that the three 
major families Methylococcaceae, Methylocystaceae, and Beijer- 
inckiaceae can be unambiguously reconstructed (Lau et al, 2013). 
Nonetheless, mxaF and 1 6S rRNA gene phylogenies differ on genus 
and species level (Lau etal., 2013). The loose coupling of mxaF 
phylogenies with 16S rRNA gene phylogenies is reflected by a low 
DNA-level similarity (about 77%) that relates to a 97% similar- 
ity cut-off on 16S rRNA gene sequence level, which is indicative 
for species. This low mxaF cut-off level even decreases when 
more species are considered (Stacheter etal, 2013) supporting 
the conclusion that horizontal gene transfer occurred between 
methylotrophs and non-methylotrophs. A frequently detected 



Table 4 | Use of amplicon pyrosequencing to analyze methylotrophic communities. 



Environment 


Gene markers 


Remarks 


Functional group 


Reference 


Soils 


Aerated soils 


mxaF, met), fael 


Amplification with adapter-less 
primers 


Bacterial methylotrophs 


Stacheter etal. (2013) 


Hydromorphic grassland soil 


pmoA 




Methanotrophs 


Shrestha etal. (2012) 


Peatland 


pmoA 




Methanotrophs 


Deng etal. (2013) 


Paddy soils 


pmoA 




Methanotrophs 


Luke and Frenzel (2011) 


Peat bog 


pmoA 


Read length > 500 nt, blended 
analysis with other genes 


Methanotrophs 


Kip etal. (2011) 


Aquatic habitats 


Lake sediments and waters 


16S rRNA genes 


Combined with DNA SIP 


Methanotrophs 


He etal. (2012) 


Water of oil sand tailings ponds 


pmoA 




Methanotrophs 


Saidi-Mehrabad etal. (2013) 


Aquifer 


Aquifer 


16S rRNA genes 


V4-V6 region of 16S rRNA 


Methylotrophs (and others) 


Lavalleur and Colwell (2013) 


Technical systems 



Methanotrophic biofilter 16S rRNA genes - Methanotrophs (and others) Kim etal. (2013) 
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mxaF genotype in temperate aerated soils is closely related to the 
methanotroph Methyloferula stellata AR4. Nonetheless, due to lack 
of congruencies between mxaF and 16S rRNA gene phylogenies 
in Beijerinckiaceae, it cannot be judged if this mxaF genotype was 
derived from a methanotroph or a non-methanotrophic methy- 
lotroph (Lau etal., 2013; Stacheter etal., 2013). Thus, a more 
in-depth analysis of bacteria that harbor MxaFI-, XoxF-, and 
Mdh2-like MDHs on the level of genomes is warranted to improve 
the understanding of the role of horizontal gene transfer and 
convergences in these organisms aiming at a more correct inter- 
pretation of mxaF, xoxF-like, and mdh2 datasets retrieved from 
amplicon-based NGS. 

FUTURE PERSPECTIVES 

In the era of metagenomics, -transcriptomics, and -proteomics, 
it is noteworthy to state that the knowledge on diversities of 
methanol-oxidizing enzymes will also facilitate the detection of 
such organisms and their metabolic pathways in "omic" datasets. 
An example is the detection of an actinobacterial MDO-like 
protein in a metaproteome of rice plants (Knief etal, 2012a). 
Since the Mdh2 or fungal MDOs have a broad substrate spec- 
trum and may utilize alternative substrates (Nakagawa et al., 2006; 



Kalyuzhnaya etal, 2008; Gvozdev etal, 2012), their role in in 
situ methanol-oxidation based solely on the detection of their 
genes is ambiguous when their activity dependent on methanol 
in situ cannot be demonstrated. Moreover, many methylotrophs 
are capable of utilization of multicarbon compounds. Thus, the 
detection of a genotype does not necessarily mean that the respec- 
tive microorganism was involved in methanol oxidation in situ. 
Hence, approaches combining gene marker-based nucleic acid 
stable isotope probing (SIP; Antony etal, 2010) or even together 
with protein SIP (Jehmlich etal., 2010) are promising to detect 
active methanol utilizers in terrestrial and other environments, and 
may allow for detection of novel oxidoreductases and microorgan- 
isms that utilize methanol in the environment. Still not resolved 
issues for a comprehensive detection of methanol utilizers are 
(a) the unresolved quantitative impact of organisms that employ 
other methanol-oxidizing enzymes than methanol-specific oxi- 
doreductases, and (b) the detection of methanol utilizers that 
dissimilate but do not assimilate methanol-derived carbon. Such 
organisms might be detectable by SIP when using their actual 
carbon substrate as a source of isotope label combined with unla- 
beled methanol and a control, in which the labeled substrate is not 
supplemented. 
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